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We study the phase transition between the high temperature Coulomb phase and the low tem- 
perature staggered crystal phase in three dimensional classical O(N) spin-ice model. Compared 
with the previously proposed CP(1) formalism on the Coulomb-crystal transition of the classical 
dimer model, our description based on constrained order parameter is more conventional, due to a 
fundamental difference between the O(N) and the dimer model. A systematic e = 4 — d and 1/N 
expansion are used to study the universality class of the phase transition, and a stable fixed point 
is found based on our calculations for large enough N. 

PACS numbers: 



The three dimensional classical dimer model (CDM), 
as the simplest model with an algebraic liquid phase (usu- 
ally called the Coulomb phase), has attracted many an- 
alytical and numerical studies [J 0, H, H| . The ensemble 
of 3d CDM is all the configurations of dimer coverings, 
which are subject to a local constraint on every site: ev- 
ery site is connected to precisely n dimers with < n < 6 
(denoted as CDM-n), and most studies are focused on 
the case with n = 1. The partition function of the CDM 
is a summation of all the allowed dimer coverings, with 
a Boltzmann weight that favors certain types of dimer 
configurations. For instance, the most standard CDM-1 
takes the following form: 



Z = E ex P(-y)> 

dimer configuration 

E = ^2 -U(n xy + n yz + n zx ). (1) 

plaqucttcs 



Tt-xy, n yz and n zx are numbers defined on each unit pla- 
quette in xy, yz and zx planes, they take value 1 when 
this plaquette contains two parallel dimers, and take 
value otherwise. When U > 0, the ground state of this 
model favors to have as many parallel dimers as possible, 
therefore when T < T c , the system develops columnar 
crystalline dimer patterns in Fig. [TJi, while when T > T c 
the system is in the Coulomb phase with power law cor- 
relation between dimer densities, which according to the 
standard Ginzburg-Landau theory can only occur at crit- 
ical points instead of stable phases. This phase diagram 
has been confirmed by a number of numerical simulations 

PBBBf. 

The nature of the transition at T c attracts most ef- 
forts. If U > 0, numerical studies confirm that this tran- 
sition is continuous, and the data suggest that at this 
transition the discrete cubic symmetry is enlarged to an 
0(3) rotation symmetry 0]. Analytically, to describe 
this locally constrained dimer system, we can introduce 
the "magnetic field" — ra^+^i, where r\i = (— 1)* is 
a staggered sign distribution on the cubic lattice. The 



number is defined on each link (Fig. [TJj), and 

n>i,fi = 1, represents the presence and absence of dimer. 
Notice that link (i, —fi) and (i—fJi, +fi) are identical. Now 
the local constraint of the dimer system can be rewrit- 
ten as a Gauss law constraint V • Bi = rji. The map- 
ping between the dimer pattern and the magnetic field 
configuration is depicted in Fig. Q] The standard way 
to solve this Gauss law constraint is to introduce vector 
potential A defined on the unit plaquettes of the cubic 
lattice, and B = V x A 0, 0|. The background stag- 
gered magnetic charge rji can be ignored in the Coulomb 
phase after coarse-graining, though it plays crucial role 
in the crystal phase. Since vector A is no longer subject 
to any constraint, it is usually assumed that at low en- 
ergy the system can be described by a local field theory 
of A, for instance the field theory of the Coulomb phase 
reads F~ J d 3 x(V x A) 2 + ■ ■ which is invariant under 
gauge transformation A — > A + V/. 

Now the Coulomb-crystal phase transition can be 
described by the condensation of the matter fields 
which couple minimally to the vector potential field A: 
cos(V9 m — A), and mathematically this matter field is 
introduced because of the discrete nature of A. Due to 
the staggered background magnetic charge rji, the matter 
fields move on a nonzero background magnetic field, the 
band structure of the matter fields have multiple minima 
in the Brillouin zone, and the transformation between 
these minima encodes the information of the lattice sym- 
metry. For instance the transition of the CDM-1 model 
is described by the CP(1) model with an enlarged SU(2) 
global symmetry 0, 0] . This field theory is highly uncon- 
ventional, in the sense that it is not formulated in terms 
of physical order parameters. It is expected that more 
general CDM-n models can also be described by similar 
Higgs transition, although the detailed lattice symmetry 
transformation for matter fields would depend on n. 

In this paper we will study an O(N) generalization of 
the CDM model. We define an 0(N) spin vector S a with 
unit length J2 a (S a ) 2 = 1 on each link of the cubic 



lattice, with a = 1 ■ ■ ■ N, and we assume that the largest 
term of the Hamiltonian imposes an ice-rule constraint 
[H for O(N) spins on six links shared by every site: 

E S l» = °- ( 2 ) 

Systems with this constraint is usually called the spin-ice. 
If N = 1, each spin can only take value ±1, therefore ev- 
ery site of the cubic lattice connects to three spins with 
+1 and three spins with —1, which is equivalent to the 
CDM-3 model. Just like the CDM, at high temperature, 
there is a Coulomb phase with algebraic correlation be- 
tween the O(N) vector S a , while the low temperature 
crystal phase is controlled by other interactions. 

In addition to the large constraint Eq. [21 we can design 
the Hamiltonian as following: 

E = E JlS i,-M S i,H + E E J 2 S Zn S i+W (3) 

Ji is a Heisenberg coupling between spins along the same 
lattice axis, Ji is a Heisenberg coupling between spins on 
two parallel links across a unit square. If J% > and 
J2 < 0, in the ground state spins are antiparallel along 
the same axis, but parallel between parallel links across 
a unit square, which is an O(N) analogue of the colum- 
nar state of the CDM in Fig. [TJi. If J\ and J2 are both 
positive, in the ground state the spins are antiparallel 
between nearest neighbor links on the same axis, as well 
as between parallel links across a unit square, which is 
an analogue of the staggered dimer configuration in Fig. 
\Tjp. Since in the J\ — J2 model Eq. [3] there is no coupling 
between different axes along different directions, for both 
cases the zero temperature ground state of model Eq. [3] 
has large degeneracy, because the spins on x, y and z 
axes are ordered independently, and the energy does not 
depend on the relative angle between S* , S? y and Sf z 
i.e. the ground state manifold has an enlarged [0(7V)] 3 
symmetry. Compared with the dimer system, the O(N) 
spin vector is not discrete, therefore we cannot use the 
Higgs mechanism as we did for the dimer model, which 
encodes the information of the discreteness of the dimer 
variables. Instead, we are going to describe the transition 
with a Ginzburg-Landau (GL) theory with order param- 
eters. The applicability of our GL theory to the original 
CDM-n will be discussed later. 

In this paper we will focus on the staggered spin or- 
der. Following the magnetic field formalism of the CDM 
mentioned before, in order to describe this system com- 
pactly, we introduce three flavors of O(N) vector field 
0^ = Sf +fl r)i with (j, = x,y, z, and now the constraint 
Eq. [2] can be rewritten concisely as 
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FIG. 1: (a), the low temperature columnar crystalline pattern 
of Eq. [T]studied in previous references; (6), the magnetic field 
distribution corresponding to this crystalline pattern (a), with 
a zero net magnetic field in the large scale; (c), the staggered 
dimer order; (d), the net nonzero magnetic field corresponding 
to the dimer crystal pattern in (c). 

Under the lattice symmetry transformation, <$> a trans- 
forms as 

T M : n -f 11 + 1, ft -> (for all fi, v), 

R^.s : fx -» -n, <j>l -> £ -> ft, (for v ^ /Li), 

: fi^v, 4>l<r+(f> a in {Iotv^h). (5) 

T M is the translation symmetry along fi direction, R^s 
is the site centered reflection symmetry, and R^ is the 
reflection along a diagonal direction. 

The mapping between Sf^ and 0° is very similar to 
the dimer case in Fig. [TJ and the staggered spin order 
corresponds to the uniform order of </>°, and all flavors 
of spin vectors are ordered. Therefore presumably 0° 
are the low energy modes close to the transition, and we 
can write down the following symmetry allowed trial field 
theory for 0° with softened unit length constraint: 

^ = E V " + * - + E 3(E V ^m) 2 + -^(6) 

When we take g — ► 00, the constraint Eq. [l]is effectively 
imposed. In Eq. when 7 = 0, the quadratic part of 
the field theory is invariant under 0(7V)xO(3) transfor- 
mation, the 0(3) symmetry is a combined flavor-space 
rotation symmetry. 7 term will break the 0(3) symme- 
try down to the cubic lattice symmetry. However, to the 
accuracy of our calculation in this paper, the RG flow 
of 7 is negligible. Therefore 7 is a constant instead of a 
scaling function in the RG equation, we will tentatively 
take 7 = for simplicity. In the limit with g — * 00, when 
r > the correlation function of 0° reads: 

m$<t>U-$) ~ rxV/.*, P^ = ^u-^. (7) 



E V ^ = °- 



(4) 



P^v is a projection matrix that projects a vector to the 
direction perpendicular to its momentum. After Fourier 
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transformation, this correlation function gives us the 1/r 3 
power-law spin correlation of the Coulomb phase. When 



r < 0, the vector $> a is ordered. 



J-4 in Eq. [5] includes all the symmetry allowed quartic 
terms of </>" : 



^ = »EEW a ] a 



^EE^iE^- ( g ) 



The it and u terms are invariant under an enlarged sym- 
metry [0(7V)] 3 , while the w term breaks this symmetry 
down to one single O(TV) symmetry plus lattice symme- 
try. As already mentioned, the ground state manifold of 
model Eq. [3] has the same enlarged [O(TV)] 3 symmetry. 
However, the w term can be induced with thermal fluc- 
tuation through order-by-disorder mechanism or we 
can simply turn on such extra bi-quadratic term energet- 
ically in the J\ — Ji model Eq. 03 Just like the J\ — J2 
model on the square lattice [9J], the quadratic coupling 
2a 0u0" with [1^1/ as well as more complicated quartic 
terms like [J2 a 0S^S]Eft 4>\4> b z\ break the reflection sym- 
metry of the system, and hence are forbidden. 

Now a systematic renormalization group (RG) equa- 
tion can be computed with four parameters it, v, w and 
r, at critical point r — with the correlation function Eq. 
[71 In our calculation we will use e = 4 — d expansion, and 
keep the accuracy to the first order e expansion. Based 
on the spirit of e expansion, all the loop integrals should 
be evaluated at d = 4, and because of the flavor-space 
coupling imposed by the constraint Eq. [2] we should gen- 
eralize our system to four dimension, and also increase 
the flavor number to /i = 1 • ■ • 4. The full coupled RG 
equation reads 
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dhd 
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2r - 6(2 + N)ur - 9Nvr - 9wr. 



(9) 



Solving this equation at r = 0, we find stable fixed point 
for large enough TV. Expanded to the order of e/7V 2 , the 
stable fixed point is located at 
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This expansion is valid in the limit e ~ 1/N 2 <C 1, and 
the fixed point values are determined by solving the equa- 
tion Eq. i]to the order of e 2 /7V 2 ~ 1/7V 6 . Close to the 
stable fixed point, the three eigenvectors of the RG flow 
have scaling dimensions 

A i = - e + °(^)> 

4448e e , 

= - e+ 567A7 + ° ( ^ 

24950e e , 

= - £+ ^7N + °^- <"> 

A3 is the largest scaling dimension, and according to Eq. 
[IT] the critical TV is N c = 44. In addition to this stable 
fixed point, there are seven other instable fixed points 
for TV > N c . For instance, in the large- TV limit there is a 
fixed point at u* = e/(247V), u* = e/(12N) and = 0, 
which has the enlarged O(TV) x 0(3) symmetry. 

In the large- TV limit, the RG equation for w is decou- 
pled from u and v, hence four of the eight fixed points 
have = 0, and all the others have = 3/ (77V). If we 
take w = w* = 3/ (77V), the RG flow diagram for u and v 
in the large- TV limit is depicted in Fig. [2] At the stable 
fixed point Eq. 1 101 r is the only relevant perturbation, 
with scaling dimension 



e + 



158e 
77^ 



(12) 



Since at the ground state all three flavors of spin vectors 
are ordered, in the field theory v should be smaller 
than 2u, which is well consistent with the stable fixed 
point in Eq. [TTJ with negative w* . This fixed point has 
positive w* , which favors noncollinear alignment between 
spins on different axes. Therefore the transition between 
Coulomb and noncollinear staggered state has a better 
chance to be described by this fixed point. 

Notice that had we included the anisotropic velocity 
7 into account, its leading RG flow will be at order of 
e 2 /TV ~ 1/7V 5 , and the flow of 7 will contribute to the RG 
flow of it, v and w at order of e 3 /7V 2 ~ 1/7V 8 , therefore it 
is justified to take 7 a constant in our calculation. When 
7 is nonzero but small, the RG flows will only change 
quantitatively. For instance, expanded to the first order 
of 7, the scaling dimensions of the three eigenvectors of 
the RG equation at the stable fixed point become Ai = 
-e A, = -e + + 2849936e ^ A, = -e 



24950e 
5677V 



30088528e7 ^ anc j scaling dimension of r becomes [r] 
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If we take TV = 1, the v and w terms are identical. 
In this case in addition to the trivial Gaussian fixed 
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FIG. 2: The schematic RG flow diagram for equation Eq. |5]in 
the large- N limit with w = u>» = Mr. The four fixed points 
A, («.,«.) = (0,0); S,(^,^f); C,(^F,-^r); and 
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C is the stable fixed point. 



point, there is only one other fixed point at v* = 2u sr — 
e/34 with 0(3) flavor-space combined rotation symme- 
try, which is the same fixed point as the ferromagnetic 
transition with dipolar interaction 1^, 11 1. In 3d space, 
the dipolar interaction also projects a spin wave to its 
transverse direction. The dipolar fixed point is instable 
against the 0(3) to cubic symmetry breaking, therefore 
when N = 1 our first order e expansion predicts a first or- 
der transition. We already mentioned that the case with 
N = 1 is equivalent to the CDM-3. However, in our GL 
formalism, in the ordered phase, the power law spin-spin 
correlation still persists if the long range correlation is 
subtracted. For instance, the fluctuation 8<f)^ — 4>^ — {<f>^) 
is still subject to the constraint ^ V M <50 M = 0, there- 
fore although the fluctuation is gapped, it still leads to 
the 1/r 3 power-law correlation. But in CDM, the ordered 
phase only has short range connected dimer correlation 
on top of the long range order [l^] • This difference is due 
to the fact that our formalism does not encode the dis- 
creteness of the dimers. Therefore one possible scenario 
for the CDM-3 with staggered ground state is that, if 
we lower the temperature from the Coulomb phase, after 
the first order transition of cf> a , there has to be another 
"Higgs" like phase transition that destroys the power-law 
connected correlation. Or there can be one single strong 
first order transition that connects the Coulomb phase 
and staggered dimer crystal directly. 

The J\ — J r 2 model Eq. [3] is invariant under cubic sym- 
metry transformation. One can turn on various types of 
cubic symmetry breaking anisotropy to this model, and 
the energy favors one of the three flavors of vector field 
(f>^ to order. Let us assume is ordered at low tempera- 
ture. To describe this transition we can turn on an extra 
mass gap m for $J and (f>y in Eq. [6l and by taking the 
limit g — > oo, at the critical point r = the renormalized 
correlation function for <A" becomes 



In this correlation function the scaling dimension [q z ] = 
2[q x ] = 2[q y ] = 2, therefore this transition is effectively a 
z = 2 transition, and J-4 is marginal according to power- 
counting. This cubic symmetry breaking situation has 
been studied in Ref. [13] . 

In this work we studied the transition between high 
temperature Coulomb phase and the low temperature 
staggered spin ordered phase in the O(iV) spin-ice model. 
Higher order e and 1/N expansion are demanded to ob- 
tain more quantitatively accurate results. The model Eq. 
[3] can be simulated directly numerically, and our RG cal- 
culation can be tested. The columnar phase with J\ > 
and J2 < is also interesting. But since the columnar 
order does not correspond to the uniform order of 0", the 
order parameter description is more complicated. We will 
study this situation in future. 

The CDM is also considered as a simple analogue of 
the spin- ice materials such as Ho2Ti20y and Dy2Ti207 
0, El El, where the Ho 3 + and Dy 3 + magnetic mo- 
ments reside on the sites of a pyrochlore lattice, and the 
ground state of these moments satisfies the same ice-rule 
constraint as Eq. [2] The Coulomb phase of the spin 
ice materials with fractionalized "monopole" like defect 
excitation has been observed experimentally 17]. The 
formalism developed in our work is largely applicable to 
the pyrochlore lattice, while the symmetry analysis and 
the number of quartic terms are different. A complete 
symmetry analysis is demanded in order to correctly un- 
derstand the O(N) spin- ice on the pyrochlore lattice. 

The author thanks very helpful discussion with Leon 
Balents, Subir Sachdcv and T. Senthil. This work is 
sponsored by the Society of Fellows, Harvard University. 
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